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Abstract 

Cooper pairing is studied in three dimensions to determine its binding energy 
for all coupling using a general separable interfermion interaction. Also consid- 
ered are Cooper pairs (CPs) with nonzero center-of-mass momentum (CMM). 
A coupling-independent linear term in the CMM dominates the pair excitation 
energy in weak coupling and/or high fermion density, while the more familiar 
quadratic term prevails only in the extreme low-density (i.e., vacuum) limit for 
any nonzero coupling. The linear-to-quadratic crossover of the CP dispersion 
relation is analyzed numerically, and is expected to play a central role in a model 
of superconductivity (and superfluidity) simultaneously accommodating a BCS 
condensate as well as a Bose-Einstein condensate of CP bosons. 
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I. INTRODUCTION 



The large-momentum divergence in the Cooper pair (CP) problem |T| was originally elimi- 
nated by a momentum-space cutoff introduced in what is now known as the Bardeen-Cooper- 
Schrieffer (BCS) model interaction which was also successfully used by BCS in the study of 
conventional low-temperature superconductors 0j. On the other hand, the short coherence 
length of some high-T c superconductors || possibly implies a stronger interfermion interaction. 
In the limit of very strong coupling one gets well-isolated "diatomic molecules" or dimers of 
fermions, as opposed to strongly overlapping CPs in the weak-coupling limit of BCS super- 
conductivity In strong coupling these dimers can conceivably undergo Bose-Einstein 
(BE) condensation (BEC). Although there is considerable controversy over the precise pair- 
ing dynamics in so-called "exotic" [0] superconductors, tracking this many-body problem from 
strong to weak coupling — known as the BCS-Bose crossover — has promoted the understanding 
of various properties of exotic materials PH^-PHP^T] . 

Here we focus on how Cooper pairing itself evolves from weak to strong coupling. We also 
study the excitation of CPs with nonzero center-of-mass momentum (CMM), which should play 
an important role in a superconducting or superfluid transition of the many-fermion system 
simultaneously exhibiting both BCS and BE condensates as in the formulation by Friedberg 
and Lee and more generally by Tolmachev [ |T9|j . 

The BCS model interaction simulates a phonon-mediated force which is effective in explain- 
ing the isotopic effect by appropriate momentum-space cutoffs placed symmetrically on either 
side of the Fermi surface However, exotic materials do not exhibit such systematic iso- 
topic effects. Also, as coupling increases the chemical potential \i decreases in value from the 
(positive) Fermi energy Ep (its value at zero temperature and zero interaction), and even turns 
negative as the so-called Bose regime is entered f§-|6||. The Fermi surface then washes out, 
eliminating any physical motivation for the BCS model interaction with Fermi-surface depen- 
dent cutoffs. Alternatively, in the BCS-Bose crossover a renormalization procedure pl| , |2"2"|| may 
be used to handle the large-momentum divergences. This leads to a renormalized dynamical 
model expressible in terms of physical observables of the system rather than ad hoc cutoffs. 

Here we derive renormalized t-matrix and Cooper equations for a pair of fermions which 
move in vacuum and in the Fermi sea, respectively. In three dimensions (3D) Cooper binding 
is expressed in terms of the two-fermion scattering length in vacuum. For a CP with a nonzero 
CMM we define a pair excitation energy as the (positive) difference between the CP binding 
energy at zero and at a finite CMM. For high fermion density and any coupling only a linear 
term in CMM dominates fE3l,E£l the CP excitation, which was in fact mentioned as far back as 
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1964 (Ref. PU| , P- 33). At any coupling and for vanishing fermion-number density a quadratic 
term dominates which is just the kinetic energy of the composite pair. The crossover from a 
linear to a quadratic dispersion for the pair excitation energy is then illustrated via numerical 
calculations. 

The CP dispersion relation enters the BE distribution function of the boson number equation 



from which T c is extracted [25|. The linear CP dispersion relation for weak coupling leads to 
novel phase-transition properties in either a heuristic p6| or a first-principles |27|] BEC picture 
of superconductivity as BE-condensing CPs. It is common knowledge that BEC is possible 
only for dimensions d > 2 for bosons with quadratic dispersion; this limitation reappears in 
virtually all BEC schemes thus far applied to explain superconductivity [|il^ ,|H|-|T9|| . But for 
bosons with a linear dispersion relation found here in weak and medium coupling, BEC can 
now occur for all c? > 1. This should be relevant in models of superconductivity encompassing 
both BCS and BE condensates |T1JT|]. 

In Sec. II the two-body problem is formulated in vacuum for a short-range, separable 
interaction. In Sec. Ill the renormalized CP equation is derived for nonzero CMM. In Sec. IV 
the CP dispersion relation in CMM is obtained numerically. Finally, Sec. V offers discussion 
and Sec. VI conclusions. 



II. TWO-BODY PROBLEM IN VACUUM 

Consider N fermions in a box with sides of length L that interact via an S-wave short-range, 
attractive (rank-one) separable potential in 3D of the form [|17j 

V pq = -(v /L 3 )g p g q , (1) 

where Vq > is the interaction strength and g p are dimensionless form factors g = (1 + 
P 2 /Po)~ 1 ^ 2 H31' wnere t ne parameter po is the inverse range of the potential. Such an interaction 



model may mimic a wide variety of short-range effective interactions: a force mediated by 
phonons, plasmons, excitons, magnons, etc. or even a purely electronic interaction. Here 
Po — > oo implies g = 1 and corresponds to a zero-range potential. The advantage of potential 
(|I|) is that many problems then yield analytic solutions. The BCS model interaction is a special 
case of ([l|) when g p is constant in the interval Ep — %ud < h 2 p 2 /2m < Ep + hujp, and zero 
otherwise, where ujp, is the Debye frequency. More realistic potentials can be approximated by 
a rank-N separable potential . 



The Lippmann-Schwinger equation for the t-matrix with potential (|l|) between two fermions 
each of mass m in free space is 
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UE) - V pq + £ V P* E _ n 2 k 2 /m + i0 t ^ E )> 

where E is the two-particle energy. For potential ([I]) the solution of Eq. (fj) is p!f 



tpq(E) 



9 P 9q 
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v ^ E -h 2 k 2 /m + i0 



(2) 



(3) 



In the limit L — > oo, the momentum sum may be replaced by an integral 

L 3 



E 



V 



d 3 k, 



(4) 



where v is the spin degeneracy. Using the zero-energy, on- shell t-matrix is given in terms 
of the (S-wave) scattering length a by f28| 



*oo(0) 



mvL 3 ' 



Eq. 



for E = then becomes 

v mL 3 
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(6) 



since #o = 1j where the iO term in the denominator is unnecessary as the sum no longer diverges 
in the small k limit. 



III. TWO-BODY (COOPER) PROBLEM IN FERMI SEA 



The CP equation for two fermions above the Fermi sea with momentum wavevectors ki and 
k 2 is given by 



2 U2 



h z k 



m 



E K + 



h 2 K 



2-| 



4m 



Ck — — E VkqCq, 



(7) 



where k = ~(ki — k 2 ) is the relative, and K = ki + k 2 the center-of-mass, momentum wavevec- 
tors, E K = 1E F — A K the total pair energy, A K > the CP binding energy, C q its momentum- 
space wave function, and the prime on the sum implies restriction to states above the Fermi 
surface: viz., |k± K/2| > k F , where kp is the Fermi wave number. For potential ([I]), Eq. (0) 
can be solved and A# determined from 
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ft 2 fc 2 /m + - 2£ F + Ti 2 K 2 /Am 



(8) 



Although the summand in Eq. (||) is angle-independent, the restriction on the sum arising from 
the full Fermi SCcl IS db function of the relative wave vector k, and therefore angle- dependent. 
The interaction strength vq in Eq. (|8]) may then be eliminated by combining with Eq. (|6|). 
Thus, in terms of the scattering length a, the renormalized CP equation is 



£' 



Wk 2 jm + A K - 2E F + h 2 K 2 /Am 
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Y h 2 k 2 /m 



muL 1 



AtxU a 



(9) 



For an attractive interaction with = 1 the sums in Eq. (0) have large-momentum divergences. 
However, each of the sums of Eq. (^) diverges in the same fashion so that the difference is 
finite. A similar renormalization of a general scattering equation can be performed [ET],2^]. 
Also, renormalized BCS gap and number equations have been used in the BCS-Bose crossover 
problem IQ-R^-O]. We shall determine the binding energy A K with Eq. (Hp employing a as 
interaction coupling parameter instead of the potential parameter Vq. In what follows variables 
are dimensionless and expressed in terms of kp or Ep = K 2 k F /2m, viz., £ = k/kp, K = K/kp, 
Ak = Ax/Ep, a = akp, etc. 

The sums in Eq. (|9]) are transformed into integrals using Eq. (|4]) . The restriction in the first 
term of Eq. ((^) arising from the full Fermi sea leads to two expressions depending on whether 
K = K/kp is < 2 or > 2, as discussed in the Appendix, 



tt/2 



dd sin 9 



1JI 



(10) 



j; /2 d9 sin 9 

r 6o 



e 9 i 



I >\n()d() / df n , 



71 

25' 



K > 2, 



where = 1 



A K /2 - K 2 /A, f3 2 K = -a&, = K 2 sin 2 9/ A + iTcos^/2, = 

1 - X 2 sin 2 0/4 + Kcos9/2, O = arcsin(2/K) < tt/2, with 9 the angle between k and 
K. To deal with the large-momentum divergences we have introduced a finite upper limit A 
and eventually let A — > oo. The momentum-space integrals in Eqs. fllt^) and (|TT| ) are easily 
performed for a zero-range interaction, i.e., with g^ = 1. For CPs with K = 0, only the first 



equation is of concern. Since in this case £o(0) 
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«o > 0, (12) 



7T I 7T 

-/3 + 1 - A) arctan — = — , a 2 < 0. (13) 

These two equations relate the dimensionless scattering length — oo < a = k F a < +oo to the 
dimensionless CP binding energy for zero CMM Aq = A^/Ep for all coupling. Transcendental 



equation (12) can be solved for CP binding for weak coupling. In this limit the argument of 



the logarithm in Eq. (|T2|) reduces to Aq/8, and one obtains 

A — > (8/e 2 ) exp(— vr/|a|) (weak coupling), (14) 



a limit first reported by Van Hove ||29|| . One can also find Ao for strong coupling or for Ao — > oo. 
Since in this limit (3q arctan (l//3 ) ~ 1, we obtain from Eq. ([T3|) 



A — > 2/a 2 (strong coupling). (15) 
Equations (g2j) and ([l|) were solved numerically to obtain the exact functional dependence of 



Ao on 1/5 , and this is compared with asymptotic forms Eqs. ([U]) and fllq). In Fig. 1 we plot 
Ak vs 1/2 spanning weak to strong coupling. One sees how the asymptotic form given by Eq. 
(|TJ]) (short-dashed line) coincides with the exact K = result in weak coupling, whereas Eq. 
(|H)|) (long-dashed curve) is also quite accurate over a sizeable region for strong coupling. 



IV. COOPER PAIR DISPERSION CURVES 



Equations (ffD) and (11) are valid for all K > and all coupling. They can be solved 
numerically for CP binding A# for any K. Before discussing numerical results we derive 
analytically the small-CMM behavior for zero range using = 1 for weak coupling in Eq. (|10|) 
which we take both for a small but non-zero K and for K = 0, and then subtract one equation 
from the other. A small-CMM expansion of the resultant equation leads to the weak coupling 
expression 

\ime K = -hv F K + 0(K 2 ) + ..., (16) 

A <1 2 

where a positive CP excitation energy Ek = (A — A#) has been defined, and the Fermi velocity 
vp is given by Ep/kp = hvp/2. The coefficient of the linear term depends only on properties of 
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the Fermi sea and not on any parameters of the potential. In contrast, the complete excitation 
energy does depend on the coupling parameter a = k F a. 

It is this excitation energy that must enter the BE distribution function in determining 
the critical temperature in a picture of superconductivity as a BEC of CPs p"7]-|T^j26|j . The 
leading term in ([16]) is linear in the CMM, followed by a quadratic term. But it is only for 
sufficiently small fermion density, i.e., when kp or Ep — > 0, and for any nonzero coupling, that 
the quadratic term dominates, viz., 

h 2 K 2 , N 

£K 2(2mj ' (17) 

the familiar nonrelativistic kinetic energy of the composite pair of mass 2m and CMM h.K. 
As mentioned, it is this dispersion relation that has been assumed in virtually all BEC studies 
of superconductivity |l^JT7H l9|. However, recent calculations of root-mean-square radii in 
two (2D) and three (3D) dimensions in the BCS-Bose crossover scheme, when compared with 
experimental coherence lengths of several typical 2D cuprates Jl5| as well as of 3D materials |l(J . 
suggest that they are describable well within the BCS (weak-coupling) regime and away from 
the Bose (strong-coupling) one. This implies that the linear approximation to the dispersion 
relation would be relevant in these cases, and that perhaps a more general description of the 
BEC of CPs for all coupling might require the exact dispersion relation. 

In Fig. 2 we display the reduced CP excitation energy Ek/A q as a function of reduced 
CMM Kj kp for zero- and finite-range potentials. Note that the CPs break up when Ek / A = 
(A — A^)/A = 1, i.e., when Ax vanishes and turns negative. These points are marked 
by dots in Fig. 2. For zero range we solve Eqs. (|TD|) and ( |TTD for g% = 1 and for typical 
values of A /E F spanning weak to strong coupling. For finite range we display results using 
g p = (1 +p 2 /Po) _1//2 w hh po = kp (i.e., range 1/po of the order of the average interfermion 
spacing ~ 1/kp). Also shown in Fig. 2 is the quadratic approximation in K as given by Eq. 
(|T7|). We have labeled the curves by A /E F as we found that the zero-range curves are closer 
to the corresponding finite-range ones than if they are labeled by l/kpa as in Fig. 1. The linear 
approximation Eq. ([16]) is valid only in the very weak-coupling and/or high fermion density 
limit. The quadratic term dominates only at vanishing density, for any nonzero coupling. For 
finite-range interaction the crossover in Fig. 2 is characterized by an inflection point with 
positive slope while for zero-range there is no such inflection point. 
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V. DISCUSSION 



The single-CP problem treated here may appear academic at first. However, it has serious 
consequences. Our CPs are taken as "bosonic" even though they do not obey (Ref. [^Df P- 38) 
Bose commutation relations. This is because for a given K they have indefinite occupation 
number since for fixed K there are (in the thermodynamic limit) an indefinitely large number 
of allowed (relative wavenumber) k values. Hence, for any coupling — and thus any degree of 
overlap between them — CPs do in fact obey the Bose-Einstein distribution from which BEC is 
determined. There have been attempts []T8| , p~9|| to formulate the superfluid and superconducting 
transition problem in a many-fermion system by accommodating both BE and BCS condensed 
phases. In these studies, the BCS (BE) condensed phase dominates for weak (strong) coupling. 
For intermediate coupling one could have both types of condensation with a certain density of 
CP bosons (fermions) available for BE (BCS) condensation. However, the full boson dispersion 
relation should be used to calculate the BEC transition temperature. 

The linear dispersion relation of a CP should not be confused with the linear dispersion 
of Anderson-Bogoliubov-Higgs (ABH) many-body excitation phonon-like modes. Collective 
modes in a superconductor were studied since the late 1950's by several workers. A more 
recent treatment for ID, 2D and 3D is available which confirms the linear ABH form 
hvpK/ yd for d =1, 2 or 3 in the zero-coupling limit. ABH phonons (like photons or plasmons, 
etc.) cannot suffer a BEC as their number is always indefinite. The number of CPs, on the 
other hand, is fixed, say, at half the number of (pairable) fermions if all of these are imagined 
paired. 

The above crossover of linear to quadratic forms of the CP dispersion relation was also 
found in 2D [ 23L]5I| so that we include both 2D and 3D cases in the following discussion of the 
BEC transition temperature T c . The 2D case is specially interesting as T c is zero for the usual 
quadratic dispersion relation of CP bosons but nonzero for linear dispersion. 

The general BEC T c formula for bosons in any dimension d and with a general boson 
dispersion relation e%_ = C S K S , for s > 0, is given p6[ for d > s by 



/ - — 

~ k 



B 



sT{d/2)(27r) d n B is/d 
2-K d l 2 T(d/s)C(d/s) 



but vanishes for d < s. Here % is the number density of bosons of mass m#, ks the Boltzmann 
constant, and ((d/s) the Riemann Zeta function. For quadratic dispersion s = 2, C s = h 2 /2rriB 
and in 3D C(3/2) ~ 2.612 ([18]) leads to the familiar T c formula T c ~ 3.31h 2 n 2 B /3 /m B k B H and 



to the fact that T c = for all d < 2. For the linear dispersion case s = 1, and consequently 
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T c = for all d < 1 and T c > for all d > 1. The latter is precisely the range of effective 
dimensionalities for all known superconductors if one includes the quasi-lD organo-metallic 
Bechgaard salts f3^.[33[ . 



Before discussing the consequences of the T c formula ( |T8| ) with s = 1 in superconductivity, 
we stress its limitations. Firstly, this T c formula with s — 1 is derived with the linear dispersion 
relation predominant for weak to moderate coupling, while the full correct dispersion relation 
should be used in general. Secondly, in deriving this formula we have taken the full momentum 
space of CP bosons so that the momentum integrals run from to oo, whereas we have seen that 
the CPs break up above some specific momentum value so the integrals should run only from 
to the breakup Kq. Thirdly, the effect of unpaired fermions in the background is ignored. 
Nevertheless, preliminary study shows |26| that once we remove these three limitations the 
result (p~8|) for s — 1 does not change drastically. 

If one assumes that all fermions are paired into CP bosons so that the boson density ris 
is n/2 with n = N/L 3 the fermion density in the normal state, Eq. (|i~8|) with s = 1 leads 
to huge values of T c ~ 10 3 K for weak to intermediate coupling — the region of interest in 
superconductivity even though T c empirically is at most about ~ 100 K. However, the number 
of paired fermions vulnerable to BEC is strongly coupling-dependent and generally |Z7 is only 
a small fraction of all the fermions so that in Eq. ( |T8[ ) is effectively much smaller than 
n/2. Thus, a realistic T c is certainly feasible. In the extreme weak coupling limit ub — > 0, 
driving the BEC T c to zero and allowing the BCS theory to be recovered from analyses as in 



Refs. [ 178 ,|T9 1 . For higher coupling n# increases so that one accommodates both BE and BCS 
condensates. One can surmise that in a realistic theory of superconductivity BE and BCS 
condensates play their respective roles. Elaborate calculations must still be performed with the 
exact CP dispersion in order to find a more accurate T c for this many-body system. 



VI. CONCLUSIONS 

The single CP problem with non-zero CMM is studied as it evolves (or crosses over) by 
varying the interfermion short-range pair interaction from weak to strong. The CP excitation 
energy is exhibited as a function of its CMM. For weak coupling the correct excitation energy 
is a linear dispersion relation in the CMM, which changes gradually to a quadratic relation 
as coupling increases and/or density is reduced to the vacuum limit. These results will play 
a critical role in a model of superconductivity that includes both BCS and BE condensates 
18|,19]. With a quadratic dispersion the BEC T c = in 2D, from which one might infer that 



BEC is irrelevant for quasi-2D cuprate superconductivity. However, even in 2D, nonzero BEC 
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transition temperatures emerge for weak and medium coupling where the linear dispersion 
relation is found to dominate [^1|,0] , thus vindicating the relevance of BEC for such materials. 



The pioneering attempts [18|fr9j developing a model of superconductivity accommodating BE 
and BCS condensates both assumed the quadratic dispersion relation for unbreakable CPs for 
all coupling. It would be interesting to reformulate those studies using the proper dispersion 
relation with a finite breakup momentum. Lastly, although our study is based on a separable 
potential we expect our conclusions on the CP dispersion to be valid more generally, and 
applicable to superconductivity, and to neutral-atom superfluidity such as in liquid 3 He and in 
trapped Fermi gases 



VII. APPENDIX 

The restriction that both fermions lie above the Fermi sea in Eq. @ can be written as 

(k/k F ±K/2k F ) 2 - 1 = £ 2 ±£Kcos9 + K 2 /4- 1 > 0, (A. 1) 

where £ = kjkp and K = K/kp- The equality leads to two pairs of roots in £, say £1,2 = — a±6 
and £ 3)4 = a ± b, where a = (K/2) cos6>, b = \Jl — (K 2 /A) sin 2 9, and 9 is the angle between k 
and K. 

For K < 2, b > a, one root of the two pairs is positive and the other negative. Thus, Eq. 
( |A lp can be satisfied provided that £ > £1,^2, £3, £4, or specifically, if £ > £o(#) = a + b. For 
K > 2 and 9 > 9 = arcsin(2/A"), b becomes imaginary and Eq. ( [A. 1[ ) is satisfied for all £. 
Therefore, there is no restriction in the integration over £. However, for K > 2 and 9 < 6 , 
b < a the pair of roots £ ii2 are both negative while the pair £ 3j4 are both positive (with £ 3 > £ 4 ). 
Consequently, in both cases Eq. ( [A. 1| ) is satisfied only if £ is in the interval [0,^' (9) = a — b], 
and in the interval [£ o (0),oo] , respectively. Equation @ is evaluated using these restrictions 
on the £ integration. 
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IX. FIGURE CAPTIONS 



1. Dimensionless CP binding energy Ak/Ef calculated for zero-range interaction using 
Eqs. (|T0D and flljD for different CMM momentum K (K/k F = 0, 1, and 3) versus the inverse 



dimensionless scattering length 1/kp a. Short-dashed straight line ([14]) holds for weak coupling 
while long-dashed quadratic one (|i~5|) is valid for strong coupling. 

2. Reduced CP excitation energy ek/Aq versus K/Uf for different reduced couplings mea- 
sured in terms of A /Ef calculated from Eqs. (|l~0|) and (|TTD for zero-range (full curves) and 
finite-range potential with p = Icf (long-dashed). Short-dashed curves are the quadratic term 
of Eq. (]TT|) . Dots denote CMM wavenumbers Kq where the CPs break up. 
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Fig. 2 




